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Local well-posedness of compressible-incompressible 
two-phase flows with phase transitions 

Yoshihiro Shibata * 


Abstract 

This paper is concerned with the basic model for compressible and in¬ 
compressible two phase flows with phase transitions The flows are separated 
by nearly flat interface represented as a graph over the N — 1 dimensional 
Euclidean space {N > 2). The local well-posedness is proved by the 

Banach fixed point theorem based on the maximal Lp-Lq regularity theorm 
for the linearized problem. 


1 Introduction 

The three states of matter are solids, liquids and gases. The flow consisting of two 
phases, which are mixed and interacting each other, is called the two phase flow. 
To analyze the two phase flow is the important problem in the field of the fluid 
machine. For example, nowadays it is known that cavitation noise and the damage 
of hard material for tarbo-machines and ship propellers are induced by impulsive 
pressures that are caused by the collapse of cloud of bubbles in the water. In fact, 
K. Yamamoto m investigated that the water jets injected from submerged noz¬ 
zle with narrow orifice were observed by a high-speed video camera and he found 
that the several times rebound of the cloud of bubbles creates very strong pressure 
pulses which cause the cavitation noise and the damage of hard material. Moreover, 
D. Rosseinelli et al [5] showed the similation of cloud cavitation collapse by the high 
performance computer. Thus, the study of the cavitation, which is described by the 
compressible and incompressible fluid flow mathematically, has new depelopment 
in the experimental fluid mechanics and computational fluid mechanics, rather re¬ 
cently. On the other hand, the mathematical approach to two phase problem with 
liquids and bubbles is rare, even when they are described by the compressible and in¬ 
compressible viscous fluid flow with sharp interface. The author knows only results 
due to Denisova [2] and Kubo, Shibata and Soga [1]. In this paper, we start with 
the modeling of the two phase problem which can be found in the usual engineering 
text books without any mathematical proof and we prove the local well-posedness 
in the phase transition case, which has not been yet treated in any mathematical 
literature as far as the author knows. 
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2 Modeling. 

Following the Priiss idea in [^, we discuss the modeling. Let be a domain in 
the N dimensional Euclidean space K.^ {N > 2) with boundary Fq. Let r2_ be 
a subdomain of with boundary F. We assume that F = d^l- C and that 
FoPlF = 0. Set = O—n_. Let ip = ip{^, t) = {ifi{^, t), - ■ ■, t)) be a function 

defined on the closure of V, for each time variable t G (0, T), ^ = (^i,..., being 
the reference coordinate system. We assume that the map ^ —>■ ip{^, t) is one to one 
for each t G (0,T)[i|. Set {dtip){^,t) = v{x,t) with x = ip{^,t), 

n±{t) = {x = ifii.t) I i G n±}, T{t) = {x = U G r}, 


and (l{t) = n_(t) un+(t). Let np(j) be the unit outer normal to T{t) pointed from 
to n+(t) and let npo the unit outer normal to Fq. Set 


[[p]] = V- — v+ (the jump of v accross T{t)) 


for any v defined on Here and hereafter, we write Moreover, 

given v± defined on H±(t), we define v by v{x) = v±{x) for x G Let 

Hr = —divpnr be the mean curvature of T(t). During our modeling, we use the 
well-known Reynolds transport theorem: 


d 

dt 



fdx 


in{t) 


dtf dx+ [[/]]v • nr(t) diy+ /v • npo diy, 


/r(t) 


^ro 


where dv represents the surface element not only of T(t) but also of Fq. In this 
orientation, we know that 


^\r{t)\ = -f Hrv-nrdi^ ( 2 . 1 ) 

« dr(i) 

Here and in the following, we use bold small letters to denote iV-vector and 7V-vector 
valued functions and the bold capital letters to denote N x N matrix and N x N 
matrix valued functions, respectively. For v = (ui,..., vn) and w = {wi ,..., wat), 
we set <v, w>=v-w = which is the usual inner product of . 

In the following, we use the following notation: 

• p : Cl{t) R+ is the mass field, 

• u : tl{t) the velocity field, 

• TT : (l{t) R the pressure field, 

• T : Q{t) {A G GLn(M.) \ = A} the stress tensor field, 

• 0 : H(t) —>■ ]R_|_ the thermal field, 

• e : (l{t) ^ M. the internal energy, 

• q : H(t) —)• the heat flux, 

t Since this subsection is concerned with the modeling, we do not care the regularity of boundary 
and the map ip. Moreover, we do not mention any integrability of functions regorously. These are 
formulated mathematically in sections 2. 
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• 77 : J7(i) —)• R the entropy, 

where, we have set R+ = (0,oo). For the modeling, we use the following Navier- 
Stokes-Fourier system of equations: for x S fl{t) ad t > 0 


dtp + div (pu) = 0; 

dt{pu) + div (pu ® u) — div T = 0; 


( 2 . 2 ) 

(2.3) 


5i(||up + pe)+div((||up+pe)u)-div(Tu-q) =0. (2.4) 

Here, for any u = (ui,..., un), u ® u is the N x N matrix whose {i,j) component 
is UiUj, and for any w = (wi ,.. .,wi\[) and S = (S'y) the divergence forms divw 
and div S are defined by 

N N N 

div w = ^ djWj, div s = (^ djSij, •. -, djSNj)- 

i=i 

During our discussion of the jump condition on T(t) and boundary condition on Fq. 
we assume that v ^ u on r(t), but v = u on Fq. 

First, we consider the mass conservation: 


- / pdx = 0. 


dt 


lQ{t) 


(2.5) 


By (12.211 and the Reynolds transport theorem, we have 

■ / pdx= dtpdx+ [[p]]v • nr(t) + / pu ■ nr„ di^ 
JQd) Jn(t) Jr(t) JTn 


_d 

dt 


in(t) 

/ div{pu)dx+ / [[p]]v • nr(t) + / 

dn(t) dr(t) dr 

/ [[P(u - v)]] • nr(t) di'. 

JTd) 


pu • npo diy 


Thus, to obtain (03), it is sufficient to assume that 


[[p(u-v)]] •nr(t) =0 onr(t). (2.6) 

In this case, p+(u+ — v) • = p-(u_ — v) • nr(j) on r(t), so that the phase flux 

j is defined by 

J = P+(u+ - v) • nr(t) = p-(u_ - v) • nr(t). (2.7) 

• When J = 0, we have [[u]] ■ nr = 0. 

• When J 7 ^ 0 and [[p]] ^ 0, we have 


[[u]] ■ nr(t) 

[[1/p]] 


( 2 . 8 ) 


• When J 0 and [[p]] = 0, j can not be decided by the velocity field u. 
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The case where j = 0 is called without phase transition and the case where j 7 ^ 0 
with phase transition. 

Next, we consider the conservation of momentum: 


dt 



= 0 . 


(2.9) 


By (12.31) and the Reynolds transport theorem, we have 


d 

dt 



I div (pu (g) u) dx + / 

Cl{t) J tl{t 

- v) - [[T]])nr(t) dv+ f 
Jr 


div T dx 


Tnro diy. 


Thus, in order that ^ pudx = 0 holds, it is sufficient to assume that 


l[[pu( 8 )(u-v)-T]]nr(t) =divrTr on r(t), 

( Tnro =0 on To 

Here, Tr is the stress tensor field on r(t). Note that div rTr = 0. We 
assume that divrTr = —(jJIrnr(t), where cr is a non-negative constant describing 
the coefficient of surface tension. 

We represent the interface condition (12.101) with the help of the phase flux j as 
follows: 


[[pu g (u - v)]] • nr(t) = p_u„(u_ - v) • nr(t) - p+u+(u+ - v) • nr(t) = j[[u]]. 
Moreover, by (123) we rewrite (ESI) as follows: 
dt(pu) + div (pu g u) = u(9tp -|- div (pu)) -|- p(dtu + u • Vu) = p(dtu -(- u • Vu). 


Summing up, we have obtained 


{ p(dtu -f u • Vu) — div T = 0 in iJ(t), 

j[[u]] - [[Tnr(t)]] =-cr7drnr(t) on r(t), (2.11) 

Tnro =0 on Tq. 

Here and in the following, for any TV- vector functions w = { wi ,..., wn ), z = 
(zi,..., zn) and scalor function /, we set w • V/ = J2j=i ^nd w • Vz is an 

Wvector function whose i th component is w • Vz^. 

Next, we consider the balance of energy. We look for a sufficient condition to 
obtain the conservation of energy: 

4(/ (f |up+pe)da;-bcr|r(t)|) = 0. (2.12) 

® Jn{t) ^ 

By E31) and the Reynolds transport theorem, we have 


d 

dt 


/ (f |u|^+pe)da: 

Jn{t) ^ 

- f [[(f |up + pe)(u - v) - (Tu - q)]] ■ nrp) diy + 
Jr(t) ^ 



q) • nro dv, 
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which, combined with (12.ip . furnishes that 


_d 
dt' 


(/ +pe)da; + cr|r(t)|) = /" (Tu - q) • npo d;/ 

Jn(t) ^ Jfo 




'm 


up + pe)(u - v) - (Tu - q)]] • nr(t) + aHrv ■ nr(t)) du. 


Thus, in order to obtain (I2.12L it is sufficient to assume that 

f [[(i|up + pe)(u - v) - (Tu - q)]] • nr(t) + cri/pv • nr(t) =0 on r(t), 
I (Tu - q) • uro = 0 on Tq. 


Since Tur^ = 0 on Tq, we assume that q • upo = 0 on Tq. By (12.71) and (12.111) . 


(2.13) 


[[||up(u-v)]] •nr(t) 

= l(|u_ - V + vp - |u+ - V + vp) = |[[|u - vp]] + j[[u]] . V 
= ^[[|u - vfl] + [[Tnr(t)]] ■ v - aHrv • nr(t) 

Since [[pe(u — v)]] • nr(t) = j[[e]], the first equation of (12.131) becomes: 

|[[|u - vp]] + j[[e]] - [[T(u - v)]] . nr(*) + [[q]] • nr(,) = 0. 

Moreover, using (lO) and (lO) . we rewrite dm) as follows: 

+ pe) + div ((^|up + pe)u) - div (Tu - q) 

= (^|up + e)dtp + p{u ■ dtu + dte) + (^|up + e)div (pu) + pu ■ (u • Vu + Ve) 

— (div T) ■ u — T : Vu + div q 

= (^|u|^ + e){dtp +div (pu)) + u • {p{dtu + u • Vu) - divT) 

+ p{dte + u • Vu) — T : Vu + divq. 

Here, we have set T : Vu = TijdiUj- Thus, we have 

p{dte + u • Ve) + div q — T : Vu = 0. 

Summing up, we have obtained 

{ p{dte + u • Ve) + div q — T : Vu = 0 in 
2 [[l“ - vP]] + j[[e]] - [[T(u - v)]] • nr(t) + [[q]] • nr(t) =0 on r(f), (2.14) 

q • nro = 0 on Tq. 

The number of interface conditions is so far not enough. To find one more 
condition, we consider the law of entropy increase: 

7 /» 

— ppdx >0. (2-15) 

dt Jh(t) 
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For this purpose, we introduce the constitutive laws in the phases. According to 
the Newton law, the stress tensor T is given by 

T = 2/iD(u) + (A — /i)div ul — ttI. 

Here, D(u) = i(^Vu + Vu) is the strain tensor field, I the N x N identity matrix, 
/i and A are the first and second viscosity coefficients satisfying the condition: 

N -2 

/i>0, A > jy (2.16) 

To prove local well-posedness, it suffices to assume that fi > 0 and A > 0. According 
to the Fourier law, the heat flux q is given by 

q=-dV6». (2.17) 


with thermal conductivity d satisfying the condition: d > 0. Moreover, the first 
law of thermodynamics tells us that the internal energy e, the entropy ry, and the 
pressure term tt have the relation: 


de = Odrj H— 


(2.18) 


We define the free energy ij) for the unit mass and the specific heat k by 


ip = e — Or], 


- ^ 
“ 56 »’ 


(2.19) 


dc dc 

respectively. We assume that k > 0. Since — = 6 and — 

or] op 

(HHD, by (1^ 


= — as follows from 


dtc + u • Ve 


dc dc 

■^{9tr] + u ■ V77) + —{dtp + u ■ Vp) 


= Oldtp + u • Vp)-div u. 

P 


( 2 . 20 ) 


In addition, 

T : Vu = 2p|D(u)p + (A — p)(divu)^ — Trdivu, (2-21) 

which, combined with the first equation of (I2.14|) . (12.201) . and (12.211) . furnishes 

p9{dtr] + u • Vp) — div {dVO) — (2p|D(u)p + (A — p)(div u)^) = 0. (2.22) 


On the other hand, we have 


5t(pp) + div (ppu) = p(5tp + div (pu)) + p(i9tp + u • Vp) = p{dtT] + u • Vp). (2.23) 
In the following, we assume that 


0 > 0, [[0]] = 0. (2.24) 

Since 9 represents the absolute temperature, d > 0 is natural assumption. While 
phase transition happens, the temperature does not change, so that [[d]] = 0 is also 
natural assumption. By (12.221) and (12.231) 

dtipp) + div (ppu) = ^{div (dVd) + 2p|D(u)|^ + (A - p)(divu)^}. (2.25) 

U 
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dt 


By the Reynolds transport theorem, (12.251) . and the divergence theorem of Gauss, 
pr] dx 

div{pr]u)dx+ j i{div (dV0) + (2^|D(u)p + (A —/r)(div u)^} dx 

+ / [[H]v • nr(t) diy+ ppu ■ npo diy 

Jr{t) Jto 

-f [[(Pd)(u-v)]] •nr(t)dz/+ / [[^V0]] ■ Ur^t) di^ + [ ^{dVO ■ nro) di' 

Jr(t) Jr(t) JTn " 


/n(t) 

-/ 

Jn{t) 


f I 

hit) 

/ ' 

Jn{t) 


d|V6i|2 

02 


/r(t) 

+ 2^|D(u)p + (A — /r)(div u)^} dx. 


Since 2^|D(u)p + (A — /r)(divu)^ > 0 as follows from (12.161) . to obtain (12.151) it is 
sufficient to assume that 


[[(Pd)(u - v) - -VO]] ■ nr(t) =0 on T{t), 
dVO ■ npo =0 on Fg. 

Moreover, by [[0]] = 0 and (12.71) . the first equation of (12.261) becomes 
j[N] - [MV0]] • nr(t) = 0 on r(t), 
which is called the Stefan law. In fact, by (12.7|) and [[0]] = 0, 


(2.26) 


(2.27) 


0 = [[{pv){u - v) - -VO]] ■ nr(t) 

= (0-0-)(u- - v) • nrp) - (p+? 7 +)(u+ - v) • nr(t) - 


(d_V0_ -d+V0+)-nr(t) 


= gbilN] - [MV0]] •nr(t)). 


Note that the Stefan law becomes the usual jump condition: [[dV0]] • nr(t) = 0 on 
r(<) provided that j = 0. 

Next, assuming that j ^ 0 and [[p]] ^ 0 and using (12.271) . we rewrite (12.141) . 
Given w, we set = w — (w • nr(t))nr(t), which is the tangential part of w 

along nr(t). Since w = (w • nr(t))nr(t) + Tnr(t)'^ and Tnr^t)'^ ' nr(t) = 0, we have 


|w|2 = |wnr(t)P + |7;r(,)W|' 
In the following, we assume that 

'T'nrto [[“ “ 


(2.28) 


(2.29) 


Especially, we have 

[[u]] = ([[u]] • nr(t))nr(t). (2.30) 

Since |wp —|vp = (w —v)-(w+v) for any w and v, by ()2.29l) (u—v)p]] = 0, 

so that by (12.71) and (12.281) 


l[[|(u-v).nr,„p]| = £||l] 


(2.31) 
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Since e = ij: + Or], we have [[e]] = [['0]] + \\Pr]\\. Thus, recalling (12.171) and using 
(12.271) and (I2.31L we rewrite the jump condition in (I2.14p as follows: 


0 = |[[|u - vp]] + j[[e]] - [[T(u - v)]] . nr(*) - [[q]] ■ nr(,) 
= j"[[^]]-[[T(u-v)]].nr(,)+j[M]. 

Moreover, noting that T± are symmetric matrices, we have 


[[T(u - v)]] • nr(t) = [[(u - v) • Tnr(t)]] 

= [[((u - v) • nr(t))nr(t) • Tnrp)]] + (u - v) • Tnrp)]]. 


By (EH), 

[[((u - v) • nr(t))nr(t) • Tnr(t)]] = j[[-nr(t) ■ Tnr(t)]]. 
On the other hand, by (12.111) and (12.291) and (12.301) . 


(u - v) • Tnr(t)]] = (u_ - v) • [[Tnr(t)]] 

= - v) ■ (j[[u]] ■ nr(t) + cri7r)nr(t) = 0. 


Thus, we have obtained 0 = j([[0]] +J^[[^]] - [[-nr(t) • Tnr(t)]]). Since j ^ 0, 


finally we arrive at the condition: 

1 „ .4 


[M] +J [[^]] - [[^nrp) • Tnr(t)]] = 0 on T{t), 


P 


(2.32) 


which is called the generalized Gibbs-Thomson law. 

Finally, we calculate Vr := v-nr(j). By (j2.7l) we have v-nr(t) = u_ •nr(t) — 

When j = 0, it follows from dZIl) and (12.291) that [[u]] = 0, so that v-nr(t) = u-nr(t). 

When J 7 ^ 0 and [[p]] ^ 0, by (12.81) we have j = 

J [[pu]] • nr(t) 

vnr(b=u--nr,)-- = ^p^. 

Summing up, we have obtained 


Vr := V • nr(t) = u • nr(t) (j = 0) 

[[pu]\ ■ nr(t) 


Vr := v • nr(t) = 


[[P]] 


(j yf 0 and [[p]] yf 0). 


(2.33) 


Next, we consider the case where j y^ 0 and [[p]] = 0. In this case, [[u]] •nr(t) = 0, 
which, combined with (I2.29L furnishes that [[u]] = 0, so that (I2.11|) becomes 


[[Tnr(t)]] = cri7rnr(t) on T{t). 


(2.34) 


To derive ()2.32|) . we assume that [[p]] yf 0, so that we reconsider the second condition 
of (|2.14p. By [[u]] = 0, [[|u - vp]] = 0. By (|2.17p and (|2.27p, j[[e]] + [[q]] • nr(t) = 

































j[['(/']]■ In addition, by (12.71) . (I2.34p . and the symmetricity of T± 


[[T(u - v)]] • nr(t) = (u_ - v) • T_nr(t) - (u+ - v) • T_nr(t) 

= (u- - v) • [[Tnr(t)]] = (u_ - v) • aHrUr^t) = 


Since j ^ 0, the second equation of (j2.14l) becomes 

M] - —Hr = 0. on T{t). (2.35) 

P- 

Noting that dtp + u • Vp = —pdivu as follows from (12.21) and recalling the 
. , de , de 

formulas: -^ = k and 7 ^ = (cf. (12.181) and (12.191) 1. we have 

Hh Hn 


^ N ,9e ^ ^ de ^ de de „ , 

p{dte + u • Ve) = p{^dt9 + -^dtp + • VO + —u • Vp) 

do dp do dp 

= pnidtO + u • V9) — —div u, 

P 

which, combined with the first eqution in (12.141) and (I2.21L furnishes that 
pK{dt9 + u • V9) — div {dV9) — (2p|D(u)p + (A — p)(divu)^) + 7r(l — -)divu = 0. 


Summing up, we have obtained the equations: for x € fl(t) and t > 0 


dtp + div (pu) = 0 , 
p{dtu + u • Vu) — div T = 0, 


(2.36) 


pnidtP + u ■ V 6 () — div (dV9) = (2p|D(u)P + (A — p)(divu)^) — 7 r(l — -)divu, 

P 


subject to the boudary condition: for a; G Fq and t > 0 : 


Tnro = 0 , dS/P ■ rirp =0 on Fq, 
and one of the following interface conditions: for a; G F(t) and t > 0 
(1) When j = 0, 


[[u]]= 0 , [[Tnrp)]] = aiFrnr(t), P] = 0 , 

[[dM9 ■ nr(t)]] =0, v • nr(t) = u ■ nr(t). 


(2) When j ^ 0 and [[p]] ^ 0, 

"Tnpt*) [[u]] = 0, j[[u]] - [[Tnr(t)]] = -CTiFrnr(t), [[9]] = 0, 

j[N] - • nr(t)]] = 0, [[^/>]] +f[[^]] - [[inr(t)Tnr(t)]] = 0, 


• nr(t) = 


[[pu]] ■ nr(t) 


[H] ’ ^ [[1/p]] ■ 

(3) When j 7 ^ 0 and p = p_ = p+ (constants). 


[[uj] • nr(t) 


(2.37) 


(2.38) 


(2.39) 


[[u]]=0, [[Tnrj] = criFrnr(t), [[ 6 *]] = 0, j[[ 6 » 77 ]] - [[dV 6 » • urp)]] = 0, 

P[[V']] - criFr = 0, V • nr(t) = u • nr(t) - 


(2.40) 
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Remark 1. Assuming that il- = and $7+ = 0, we have the one phase problem. 
In this case, as boundary condition on Fq, we have 

Tnro = criFronr(t), dV6» • npo = 0 on Fq. 


3 Problem 


The problem of this paper is concerned with the compressible and incompressible 
two phase flow separated by a nearly flat interface with phase transition. Let ho{x') 
be a given function with respect to x' = (a;i,.. ., xn-i) and we set 

= {a; = (cci,..., xn) G | ±(a;Ar — ho{x')) > 0 for x' G 
F = {a: G \ Xn = h{x') for x' G 


In this case, 17 = and Fq = 0. Let h{x', t) be a unknown function and we assume 
that the time evolutions of domains fl± and the surface F are given by 


^±(0 = = (a^i) ■ • ■)G I ±(a;iv — h(a;', t)) > 0 for x'G R^ ^}, 

F(t) = {cc G R^ I Xn = h{x',t) for x' G R^“^}. 


In this case, nr(t) = (—V'/i, 1)/sjl + |V'hp, with V'h = {dih,... ,dN-ih) {dj = 
d/dxj). Moreover, (p{x,t) = {x',xn + h{x',t)), so that v • nr(t) = dt'^ ■ nr(t) = 

dth/,/Y+WW- 

In view of (I2.2IL (I2.36p . (12.371) and (12.391) . our problem is given as follows: 


For X G f7+(t), t > 0, 

f p+{dtM+ + u+ • Vu+) - DivT+ = 0, dtp+ + div(p+u+) = 0, 

I p+K+(9t0+ + u+ • V0+) — div ((i+V0+) — T+ : Vu+ — —divu+ = 0 

I P 

and, for a: G f7_(t), t > 0, 

f p*_(9tu_ + u_ • Vu_) — DivT_ = 0, divu_=0, 

1 p^-K-idtO- + u_ ■ V6»_) - div (d_V6l_) - T_ : Vu_ = 0 


subject to the jump conditions: for x G F(t) and t > 0, 


[[-]]j^nr - [[Tnr]] = -crffrnr, 

P 

II 

o 


Am-[{d^]] = o, 

P] = o, 

(3.3) 

[M] + [[^]]j^-[[inr-Tnr]]=0, 

n . [[H] • (-V'h, I) 

IHl ' 


[[pu]] • nr(t) 

[M] 


and the initial conditions: 

(p+) u+, 0+)|t=o = (p*+ + po+) Uo+, 0* + 00 +) in fi+, 
(u_,0_)|t=o = (uo+,0* + 00+) in D_, hlt=o = ho on F. 
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Here, p*±, 0* and cr are positive constants describing the reference mass densities 
of the reference temperature of il± and the coefficient of the surface tension, 
respectively. Moreover, T± = S± — 7r±l with 


S+ = S+(u+,p+,6»+)) = ^+D(u+) + (A+ -/r+)divul, 
S_ = S_(u_,0_) = ^_D(u_). 


Here, d+ = d+{p,6), /r+ = ^+(p, 0), A+ = A+(p, 0), k+ = K+(p, 0) are positive C°° 
functions with respect to (p, 9) G (0, oo) x (0, oo), and p) and ? 7 +( 0 , p) are real 

valued C°° functions with respect to {p,9) G (0,oo) x (0, oo), while d_ = d-{6), 
P- = P-{9), K- = K-{6) are positive C°“ functions with respect to 0 G (0, oo), 
and ip-{9) and ri-{9) are real valued C°° functions with respect to 0 G (0, oo). 
Moreover, we also assume that 7r_|_ is given by 7r_|_ = P+{p,0), where P+ is some 


C°° function with respect to (p, 9) G (0, oo) x (0, oo) such that 


dP+ 


> 0 for any 


ip,9) G (0,oo) X (0,oo). 

The main purpose of this paper is to show the local wellposedness of problem 
(EH), (EH) and (El in the maximal Lp-Lq regularity class under the assumption 
that p*± and 0* satisfy the condition: 


P*- P*+t 


t{j-{9^) - V’+(p*+, 




(3.5) 


To state our main result, we transform T{t) to the flat interface. Set 

= {x = {xi,... ,xn) G I ±X]S[ > 0}, Mq = {x G I xn = 0}. 


We transfer the problem given in domains H±(t) to that in U R(Y with 

interface . Let h{x', t) be a function appearing in the definition of r(t) in (13.11) . 
Let H{x,t) be a solution to the equations: (1 — A)H = 0 in R^ with H\xj^=o = 
h{x',t), where AH = To prove the local well-posedness, we assume 

that ho is small enough, so that we may assume that 

d 1 

1 + —— H{x,t) > - for any x G R± and t G (0,T). (3.6) 

oX]\i 2 

If we consider the transformation: 


VN = XN + H{x,t), yj = Xj {j = I,. .. ,N -1), (3.7) 

then by (13.611 and r{t) are transformed to R^ and R^, respectively, because 

VN = h{y',t) when xat = 0 and = 1 + > i. 

Let u±, P+, 7r_ and 9± satisfy problem (13)21) . (13.311 and (13.411 . Set 

u±{x,t) = u±{x',XN + H{x,t),t), p+{x,t) = p+{x',Xn + H{x,t),t), 

TT_{x,t) = TT-{x' ,xn + H{x,t),t) - P+(p*+,6(*), 

9±{x,t) = 9±{x ',xn + H{x,t),t), p*+= p+(p*+,6(*), A*+= A+(p*+,6(*), 

K*+= «:+(p*+, 0*), (i*+= d+(p*+, 0*), p*_=p_(6(*), = k_(6(*), 

d*_ = d_(6»*), /i+= p+(p+,0+), A+= A+(p+,0+), p-= p-{9-), 
k+= K+{p+,9+), k-=k-{9_), d+= d+{p+,9+), d-= d-{9_), 

P+ = P+ — P*+, A± = A± ~ P*±: “ ^*+1 k± = k± — K*±, 

d+ — d+ d^^ . 
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Setting Hq = dtH, Hj = djH (j = 1 ,..., TV), we have 


Hn 


{dtf){x',XN +H{x,t),t) = dtf{x,t) - dNf{x,t) 

i + Un 


{djf){x',XN + H{x,t),t) = djf{x,t) - 


H, 


l + H 


N 


-On fix,t) {j = l,...,N). 


( 3 . 8 ) 


In the following, we set 

H, 


K, = 


1 + Hn 
By ( 13.811 we have 


(j=0,l,...,iV), K = (iVi,...,iVw), Ko = (iVo,K). 


Vtt- = QVtt- = 


/I 

0 •• 

• 0 

K, \ 


f dx^- \ 

0 

0 

• 1 

Kn -1 



VO 

0 •• 

• 0 

1+Hn / 




and Q ^ is given by 




0 

• 0 

-i?l \ 



/O • 

• 0 

-Hx \ 

Q-^ = 

0 

0 

• 1 

—Hm-1 

= I+Qi 

with Qi = 

0 • 

• 0 

—Hn-x 


VO 

0 

• 0 

1 + Hn) 



VO • 

■ 0 

Hn ) 


By ( 13.811 we have 


div u± = div u± + Vdiv 
1 


1 + Hm 
1 

1 + Hm 


{divu± - /_(u±,iV)} 
{divu± — divf„(u±,iV)} 


( 3 . 9 ) 


with 


N N-1 

ydiv(u±,I?) = -'^KjdNU±j, f-iu±,H) = '^iHNdjU±j - HjdNU±j), 

T=i i=i 

Af-l 

f_(u±, iJ) = —{Hnu±i, ..., Hnu±n-i, — ^ HjU±j). 

i=i 

For any N x N matrix of functions G = ^(gi,..., g^v), where denotes the 
transposed M, with ^g^ = {gn, ..., giN), by (13.911 we have 

DivG = DivG + VDiv(G,iV) (3.10) 

with Voiv (G,iV) = "^(^div iii,H), t4iv iiN,H)). Moreover, we set 

D,j{u±) = D,j{u±) + VDijiu±,H), D(u±) = D(u±) + VD(u±,iJ), (3.11) 
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where Voi^ (u±,i?) = —{KidNU±j + KjdNU±i) and Yr>{u±,H) is the NxN matrix 
whose (z,j) component is Vd-^ (u±, if). 

Under these preparations, we see easily that problem (13.21) . (13.31) and (13.41) is 
transformed to the following problem: 


r dtp+ + v+ • Vp+ + / 5 +(divu+ + Udiv {u+,H)) = 0 
J p^+dtu+ - Div S*+(u+) = F+ 

[ p^+K^,+dt0+ - (i*+A0+ = Fe+ 

{ Pf-dtii- — Div S*_ (u_) + Vif- = F_ 
div u_ = /_ = div f_ 

Pf-Kt^-dtO- — d*_A0_ = Fe- 

' AAr(U-)|- - p*+DiN{u+)\+ = G, 

T_-T+ = aA'H + Gn, P^IT. - = Gn+i 

^ U-i\- - u+t\+ = Ki 

0-\- - ^+1+ = 0. d^-dN0-\- - d^+dN0\+ = Gg 

dtH - - u_j,\_ -= Gh 

. ^ P*— P*+ P*— P*+ ' 

j (P+, u+, ^+)|t=o = (po+, uo+, 00+) inR)^, 

I (u_, 0 _)|t=o = (uo-,00-) H\t^o = Ho on 

where i = l,...,A — 1, and we have set 


in R^ X (0,r), 


in R^ X (0,r), 


on R^ X (0,r), 


(3.12) 


= {p^_Dnn{u-) - (3.13) 
T+ = {p*+DNNi'i^+) + (A*+ — /r,+)div u+), 

N 

V+ = (m+1, • ■ • , U+N-1,U+N — Kq — ^ KjU+j), 

1=1 

S*+(u) = ^*+D(u) + (A*+ — /r,+)divul, S*_(u) = /r,_D(u), 

Po+{x) = po+{x',XN + Ho{x)), uo±(x) = Uo±(a;',a;Ar + Ho{x)), 

0o±{x) = eo±{x',XN + Ho{x)). (3.14) 

Here, /|±(a;o) = lim f{x) for Xq £ R^. Moreover, Hq is a function satisfying 

the equations (1 — A)Hq = 0 in R^ and 1 ?oUn=o = do, and the right-hand sides in 
(13.121) are defined by the following formulas: 

F+ = F+{p+,u+,H) = 

- p+dtu+ + p+{KodNii+ - u+ • Vu+ -|- (u+ • K)clAru+} 

+ Div (m+D(u+) + A+Vd(u+, iJ)) -h Vd(A+(D(u+) + Vn{u+,H))) 

4-V{(A+ -/x+)divu+ -4 (A+ -/t+)Udiv {u+,H)} 

+ VDiv((A+ - /(+)(divu+ -4 Udiv (u+,i7))I) - QVP+(p+,0+), 

N 

F0+ = F0+{p+,u+, e+,H) = -d.+ Y, dj{K,dNh) 

1=1 
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N N 

+ Y, djid+{dj0+ - K.dNh)) + Y KjdN{d+{dj9+ - KjdN^)) 

j=l 3=1 

- {p+k+ - p^+K^,+)dtO+ + p+k+{KodN0+ - u+ • V0+ + (u+ • K)dNO+) 

+ 2/i+|D(u+) + Vd(u+, i?)P + (A+ - A+)(div u+ + (u+, H))^ 

+ P+(/5+,0+)(l - -::^)(divU 4 . + ydiv (u+,-ff)), 

P+ 

F_ = F_(u_,iJ) = -Qi(p*_i9tu_ - /i*_DivD(u_)) 

- (/ + Qi){p,_(iiro'9AfU_ - U_ • Vu_ + (u_ • K)i9ArU_) 

+ Div (m-D(u_) + A-VD(u_,i/)) + VDiv (A-(DivD(u_) + Vd(u_, H)))}, 

N-l 

f=f-{u-,H) = Y {HnOjU-j - HjdNU-j), 


3=1 


N-l 


f_ = f_(u_,i7) = -{HnU- 1 , .. .,HnU-n-i, - Y 


3=1 


N 


= Fe_(u_,0_,i/) = -p^^k-dtO- - d^-Y^ji^j^Nd-) 


3=1 


N 


+ p^-k-iKodNO- - u_ • V0_. + (u_ • K)dNO^) + Y - K^BnO-)) 


3=1 


N 


+ Y^i^N{d-{d,9_ - K^dNO-)) + 2fi_\D{vi^) + Yn{vi-,H)\^, 

3=1 

Gi = C?i(/5+, u±, H) = 

- {(/i_Aiv(U-) + A-VDi^(U-,i/))|- - (/i+AAr(u+) + A+VD,^(u+,i7))| + } 
N-l 

+ Y^^^^^^N{.Dij{n-) + VDii{-»-,H))\_ - fi+{D,j{u+) + VD,^{u+,H))\+} 


3=1 


N-l 


- i^iH)[Y^^3^'>^d-iD^Jiu_) + Vd,^{u-))\- -/i+(Aj(u+) + VDy(u+))|+}, 

3=1 

- {p-{Dnn{u-) + Vd^^(u_, i/))|_ - p+{DNjiu+) + Vd^^(u+, i/))|+}], 

Gn = GNip+,ii±, H) = 

- iJl-DNNiu-) + P-yDMNi<^-,H))\- - {fl+DNNiu+) + P+VDMNii^+,H))\ + 

+ {(A+ — /i+)divu+ + (A+ — /t+)Fdiv (u+, -ff) — (P+(p+, 0+) — P+(p*+, 6i,+ ))}| + 

N-l 

+ Y i93H){P-{DNj{u-) + Vd^j{u-,H))\_ - fi+{DNj{u+) + Vd^.(u+, P))|+) 


3 = 1 


'{(> 


N-l 


yTTWH\ 

1 1 


1 _ 

' (1 + I V'P|2)3/2 


} 


+ (-(«-Ar|--«+Ar| + f(l + |V'iJn, 

Vp*_ P+I+/ 
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Gn+1 — GAr+i(p+, U±, -ff) — 

-(-f+(P+)^+) “ P+{P*+7(^*+)) — (^r-j-)-f’+(P+)^+) 

P*- ^P+l+ P*+^ 

+ ('0-(^-) - ■0-(6'*))l- - (V’+(/3+,0+) - V'+(P*+,^'*))I+ 

+ (V>-(0-)|--V'+(p+,0+)| + )|V'i7|2 


+ - 1" - I ) (-- I ) (w-wl- - M+Jv| + )^(1 + |V'-ffp)^ 

2Vp*_ p+\+J\p^_ P+I+/ 

- i-/i_DjVAr(u-)|--(/i+£>iVAr(u+) + (A+ - /2+)divu+)|4 .1 

^ P^— P*-\- ^ 

+ -) {A+(^AfAf(u+) + ^£>jviv(u+,^^)) 


+ (A+ - A+)(divu+ + Vdiv iu+,H))}\ + 


N-l 

A ._1 ^ P^- 


- ^{D,,iu+)+VD,,in+,H))\ + }id,H)id,H) 
P-\- ^ 

N-l . 

+ 2 E I — {n~) + VD,An-,H))\- 

i=i P*- 


- ^(AAr(u+) + VDi!^iu+,H))\ + \{d^H) 

p+ J 

_ ,v.,^,2r-/ {d,H){d,d,H) ^ 

' ' L l^i + |V'i/|2 Z^^{1 + \V'H\^NP 

+ (—- -Pj-y\u-N\- - u+m\ + N1 + \V'H\^) 

Vp^_ p+l+z 

A^-1 

+ Yyp-{D^N{u-) + VDi^iu_,H))\_ - fl+{D,N{u+) + VD,r,{u+,H))\ + }i^^H) 

- {/i_DArAr(u_) + yD„„(u+,i?))|- -/i+(DArAr(u+) + yD„„(u+,iJ))| + }, 

Ki = K,{u±,H) = -{diH){u_N\- - M+iv| + ), 

Ge = Gg{p+, u±, §±, H) = 

{i + \v'H\^){u.M\--u+N\+)(—--:y-y\e-p-{e-)\--hp+{p+N)\+) 

Vp*_ p+\+J 

- {d-{V9- - 'KSnO-N - d+{V9+ - 'KdNO+)\+) ' {-VH,l) 

+ ((i*_V'0_|_ - d^+V9+\+) ■ V'H + (fi*_ajv0_|_ - d*+5Af0+|+)K • (-V'i7,1), 
Gh = Gh{p+,u±, H) = 

( --) (p*-U_Ar|_ — /9*+M+Jv| + ) + - ^U+n\ + 

Vp*_— /9+I+ p*_—p*+/ Pt.- — p+ 


P*- - P+\ + 


YidjH)u_,\. + ^±^Yi9jH)u+,\ + . 
J=i P*- ~ P+ J=i 
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The phase flux j is eliminated by using the formula: 

j = (u_jv|- - u+7v|+)(— - -)'Vl + |V'i/P on < X (0,r). 

\p^- p+\++p^+J 

Moreover, we have used the formulas: for a; S r(t) and t > 0 
7r_ - P+ + (A+ - /r+)div u+ 

N-l 

= -aHr- [[-]]j^ - ^ [[pDMj]]idjH) + [[pDnn]], 

^ i=i 

Prnrp) = {div'( ^ }(-V'g, l)/x/l + |V'gp 

where V'iJ = [diH ,..., Bn-iH) and div'v' = djVj for v' = (ui,..., u^v-i)- 

Since T- —T+ = aA'H + Gn and p^lT- — p~^T+ = Gn+i are equivalent to 

r±= A'g+ (p-^Gn-Gn+i), (3.15) 

P^— P=(=+ P^— P*-\- 

the compatibility condition for problem (15.121) is 
divu_o = f-{u-o,Ho) = divf_(u_o,Po) in 

P*-D,n{u-o)\- - Ai*+Aw(u+o)|+ = G'*(po+,uo±,^lo) (* = 1,. ■., fV - 1), 

uo-i\- - Mo+i|+ = Ki{u±,Ho) {i = 1,... ,N - 1), 

So- 1- — ^0+1+ = 0, 

d^-dN0o-\- — d*+clAr0o+|+ = Ge(po+) uo±, 0o±) ^lo)j 

(/r*+ZlArAr(uo+) + (A*+ — /r*_|_)div Uo+)|+ =- — - A'h^ 

P*— P*-\- 

H—-— ^ {p^1Gn{p+, uo±, i?o) “ GN+iip+, uo±j So±,Hq). (3.16) 

P*— P*+ 

The following theorem is the main resnlt of this paper concerning the local well- 
posedness of problem (I3.12p . 

Theorem 3.1. Let 1 < p,q < oo with 2/p + N/q < 1. Assume that and d* 
satisfy the condition (IT5]) . Then, given any positive time T, there exists an e > 0 
such that problem (13.121) admits unique solutions p+, u± and 9± with 

p+ € Wpi0,T),L,iRf))nL,{{0,T),W^iR^)), 
{u±,e±)ewp{o,T),Lg{m^))nLpi{o,T),w^{R^)), 

H € Wp{0, T), W^{R^)) n Lp((0, T), Wj^iR^)) 
for any initial data 

Po+€W^{Rl), (iio±,0o±)ei?2a-i/P)(K^), Hoew^-p/^iR^) 
satisfying the smallness condition: 

l|P0+|lwi(R(/) + X! + II^o|Ib3-i/p(jj;v) < e 

f.=± 

and compatibility condition (j3.16ll . 
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Here and in the following, Lq(IR^) and denote the usual Lebesgue 

space and Sobolev space of order m in the Lq{M.^) sense, while || • ||j;,^(kN) and 
II ■ II denote their norms, repsectively. For the Banach space X, Lp((0, T), X) 

and Wp{{0,T),X) denote the Lebesgue space and the Sobolev space with values 
in X, while || • ||Lp((o,T),jf) and || ■ || w™((o,T),jf) denote their norms, respectively. 

) denotes the real interpolation space (bF“(]R^), (R^))6i,p with 

real interpolation functor (•, •)e,p and 0 < 0 < 1. 

Remark 2. (1) The mathematical study of the compressible and incompressible 
two phase problem is quite rare as far as the author knows. First Denisova [2] 
studied the evolution of the compressible and incompressible two phase flow with 
sharp interface without phase transition under some restriction on the viscosity 
coefficients. Recently, Kubo, Shibata and Soga [1] studied the same problem as 
in [2] without any restriction on viscosity coefficients in case of without surface 
tension and without phase transition. This paper is the first manuscript to treat 
the compressible and incompressible two phase problem with phase transition 0. 
The incompressible and incompressible two phase problem with phase transition 
was studied by J. Priiss, et al. [aiziiH]. 


4 Maximal Lp-Lq regularity 

In the following, we assume that N < q < oo in view of the Sobolev imbedding 
theorem: IkllLoo(n) < C'||u||wi(o) with D = and D = R^. To solve problem 
(I3.12p . we use the maximal Lp-Lq regualrity for the parabolic equations. From this 
point of view, we represent pj^ by the integration along the characteristic curve 
generated by vJH to eliminate from the first equation of (I3.12L which is the 
hyperbolic equation for 

Given function / defined on the Lions extension Ext[/] of / is defined by 


Ext[/](a::,t) 


3/(F, -XN, t) - 2/(x', -2xn, t) 


for xat > 0, 
for xat < 0, 


Set w+ = Ext[u+] and v = (w+i,..., w)+jv-i, w+tv — Kq — Note 

that V = v+ on R^. We assume that 


T 

l|Vv(-,f)||i_^(RiV) dt < ei (4.1) 

with some small positive constant ei > 0. We use the usual fixed point argument 
to solve the nonlinear problem and in this argument we keep the situation where 
u_|_ and H satisfy (14.ip . 

‘Modeling and the main results of this paper were announced in the abstract of 39*^ Sapporo 
symposium on PDE at Hokkaido University (cf. El). 

tTani m represented the mass density with the help of the velocity field to prove the local 
well-posedness of the Navier-Stokes equations describing the compressible viscous fluid flow (cf. 
also |16l I18| ) It was also suggested by J. Priiss to the author to represnt p+ by and H with 
the help of the equation of balance of mass when the author visited Halle university in the early 
of April, 2014. 
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Let ^ be the solution to the Cauchy problem: 

|(?7,0) = ?? e 

According to Strohmer m, we choose an ei > 0 so small that the map: 77 i—^ is 
bijective on for any t S [0,T]. We denote its inverse map by f) = We 

look for p_|_ satisfying the equation: 

9t/5++ V • V; 0 ++/5+(divw++ Vdiv (w+,iJ)) = 0 inR^x(0,T). (4.2) 


Since 

= (atp+ + V • Vp+){i{r],t),t) = 
with g = -(divw+ + Cdiv (w+,iL)), defining p+{^,t) by 

p+(C,t) = {p,+ + po+{p))e- (4 3) 


with 77 = where po+{g) = Ext[/3o+] to we see that p+ is a required 

function satisfying (l4l^ with p+(^,0) = p*+ + po+iO in ®+- 

Inserting the formula of given in (H31) into the right-hand sides: F+ = 
F+ (p+, u±, iJ), A,+ = F,+ (p+, u±, 0+, iL), G, = G, (p+, u±, iL) (j = 1,..., IV + 1) 
and Gg = Ge(/5+, u±, iL) in (I3.12I1 . we have the interface problem of the final 
form, which is a quasilinear parabolic equation. As the linearized problem, we have 
the decoupled two systems. One is the Stokes equation with interface condition: 


p*+5tu+ -DivS»+(u+) = f+ 
p*_cltu_ — DivS*_(u_) -f V7r_ = f_ 
div u_ = /div = div fdiv 


in X (0, T) 

in X (0, T) 


(4.4) 


subject to the interface condition: for x S Rq and t S (0, T) 


/r*_AAr(u-)|_ - /i*+AAf(u+)|+ = g^ (i = 1,... ,fV - 1), 

{p^-Dnn{u-) -7r_)|_ = alls’ll + g^ 

(p.*+Z?ArAf(u+) + (A*+ — p,,+)divu_|.)|_|_ = a+A'H + gN+i, (4 3 ) 

77—i|— TX-i-fl-i- — hi {% — 1, . . . , M 1), 

dtH - ( -- u-N -- u+n) = d 

and the initial condition: 

u±|t=o=uo± in R±, H\t=o = Ho inR^, (4.6) 

where we have set a± = p±a{p^:- — p,+ )“^ and we have used the equivalent relations 
(13.151) . Another is the heat equations with interface condition: 


p*+K*+i9t0+- (i*+A6(+=/+ in]R(^x(0,T) 

— d*_A0_ = /_ in x (0,T) 
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subject to the interface condition: for x S Rq and t S (0, T) 


6»_|_ - 6»+|+= 0, d^+dN0-\-- d^.+dN0+\+= g (4.8) 

and the initial condition: 

0±\t=o = 0o± onR±. (4.9) 

We have the followins: theorem about the maximal L„-L„ reeularitv for problem 

gai, US), dMi. 

Theorem 4.1. Let 1 < p,q < oo and 0 <T < oo. Assume that ^ p*+- Then, 
for any initial data Uo± G ^^^^(R^) and Hg G i3g_p^^^(R'^), and right-hand 

sides of (031) and (031) 

f± G Lp((0,T),L,(R^)), /div G Lp{{0,T),W^{Rff)), fdw G W^{{0,T), Lg{R^)) 
d G Lp((0,T), W2(R^)), g, G Lp((0, T), W,1(R'^)) n ^^((O, T), ^-^(R'^)) 
h, € Lp{{0,T),W^{R^))nW^{{0,T),Lp{R^)) 

for i = 1,..., TV + 1 and j = 1,..., N — 1, satisfying the compatibility conditions: 

div Uo_ = /_ |t=o = div fdw |t=o in K-, 

^*_AAr(uo-)|-- M*+AAf(uo+)|+= gj|t=o (i = 1,...,A^-1) on R^, 
(/r*+ZlArAr(uo+) + (A*+ — /r*+)divuo+)|+ 

= cr+A'iVo + 5Af+i |t=o on Rq , 

uo-t\-- uo+^\+= h^\t=o {i = l,...,N-l) onR^. 

then, problem 031), 031), 011) admits unique solutions u± and H with 

u± G Lp{{0,T),W^{R^)) n W^{{0,T),Lg{R^)), 

H G Lp((0, T), W3(R^)) n Wpi((0, T), W^{R^)) 

possessing the estimates: 


l=± 

+ ll'9t-ff||Lp((0,t),W'2(RiV)) + ||if||Lp((0,t),W3(RW)) 

<Ce'*^{J2i\\noi\\^.a- ) + l|f^llip((0,t),i5(Rf))) + ll/div llLp((0,t)W5HR") 

i=± 

AT+l 

+ l|fdiv IIlp((0,T),L,(R^)) + X! (II^*II^p(( 0-‘).W",HR")) + ll^*5i|lLp((o,t).M/,-bR'^))) 

i=l 

N-1 

+ E (ll^ilUp((o,t)W|(R'^)) + ll'^t^jllip((o,t),i5(R''))) + IMIlLp((o,t)W|(R"))} 

i=i 

for any t G (0,T) with some positive constants C and 7 independent oft and T. 

And also, we have the following theorem about the maximal Lp-Lq regularity 
for problem (14.71) . (14.81) . (14.9|) . 
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Theorem 4.2. Let 1 < p, q < oo and 0 < T < oo. Then, for any initial data 
9 q± £ Bq^p and right-hand sides 

/± Gip((0,r),L,(R^)), gGLp((0,r),TT,i(K^))nTTpi((0,r),TT-i(R^)) 

satisfying the compatibility condition: 

[[6»o]] = 0, d^-dN0o-\-- d*-dN0o+\+= g\t=o onM.^, 

problem (lO) and (03 admits unique solutions 9± with 

9± G Lp{{0, T), )) n Wpi(( 0 , T), L,(M^)) 

satisfying the estimate: 

t=± 

< C'^*{'^{\\9oe\\B^^a-i/p)(^N) + ll/dlLp((o.t),i,(Rf))) 

i=± 

+ ll5llLp((o,t),iyi(R'^)) + ll'5‘5llLp((o,t).M^-i(RN)))} 

for any t G (0,T) with some positive constants C and 7 independent oft and T. 

Remark 3. (1) The proof of Theorem 14.Il ls given in [T3]. The proof of Theorem 
I4.2l is found in [3], but it can be proved by using the same argument as in the proof 
of Theorem 03] in |14) . 

(2) Theorem 13.11 is proved with the help of Theorem 14.11 and Theorem 14.21 the 
Banach fixed point argument and, some bootstrap arguments. The argument is 
quite standard, so that we may omit the proof of Theorem 13.II fcf. Priiss [S]). 

5 7?.-bounded solution operators 

To prove Theorem 14.11 we consider the following generalized resolvent problem: 

p*+Au+— Div S*+(u+) = f+ in 

p*_Au_ — Div S*_(u_) + V 7 r_ = f_, div u_ =/div = div fdiv in (5.1) 

subject to the interface condition: for x £ 

/r,_AAr(u-)|_ - /i*+AAf(u+)|+ = g^ (i = 1,... ,1V - 1), 
(/r,_DArAr(u_) - 7 r_)|_ = a-A'H gN, 

(/r*+ZlArAr(ud_) + (A*+ — /r,+)divu_|-)|+ = a+A'H + gN+i, 2 ) 

— hi {i — 1 ,..., M 1), 

\ TT ( P^— P*+ \ 7 

XH — I - U-N - = Ct, 

\ Pufi— P^-\- P^— P*-\- ^ 

which is corresponding to the time dependent problem (14.4L (14.51) . (14.61) . 

Before stating the main result of this section, we first introduce the definition 
of 7?.-boundedness and the operator valued Fourier multiplier theorem due to Weis 

m- 
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Definition 5.1. Let X and Y be two Banach spaces with norms || ■ ||x and || • Hy, 
respectively. A family of operators T C £(X,Y) is called 72.-bounded on C{X,Y), 
if there exist constants C > 0 and p G [1, oo) such that for any n G N, C T, 

{/j}j=i C X and sequences of independent, symmetric, {—1, l}-valued 

random variables on [0,1] there holds the inequality: 

a l n 1 pi ^ i 

- j=i “'0 j=l 

The smallest such C is called 7?.-bound of T, which is denoted by 'R-c(x,y){'T). Here, 
£(X, Y) denotes the set of all bounded linear operators from X into Y. 

Let X>(]R, A) and 5(K., A) be the set of all A valued (7°“ functions having com¬ 
pact supports and the Schwartz space of rapidly decreasing A valued functions, 
respectively, while 5'(R, A) = £(5(IR, C), A). Given M G Li locfK \ {0}, A), we 
define the operator Tm : A) ^ 5'(R, Y) by 

= (A[</)] gP(R,A)), (5.3) 

The following theorem is obtained by Weis m- 

Theorem 5.2. Let X and Y be two UMD Banach spaces and 1 < p < oo. Let M 
be a function in C'^(K \ {0}, C{X, A)) such that 

^c(x,y)({(T^)^A^('r) I T G K \ {0}) < k < oo (£ = 0,1) 

with some constant n. Then, the operator Tm defined in (lESD is extended to a 
bounded linear operator from Lp(K,X) into Lp(R, A). Moreover, denoting this ex¬ 
tension by Tm, we have 


II'7m||£(Lp(r.x).Lp(r,v)) < Cn 

for some positive constant C depending on p, X and A. 

Remark 4. For the definition of UMD space, we refer to a book due to Amann [T]. 
For 1 < q < oo, Lebesgue space Lq{TL) and Sobolev space are both UMD 

spaces. 

Theorem 5.3. Let 1 < q < oo and 0 < e < tt/ 2. Set £e,\o = {A G Ee | |A| > Aq} 
(Ao > 0) with Ee = {A = 7 -I- G C \ {0} | | arg A| < tt — e}, and 

Xq = {(f+,f_,/div,fdiv,g,h,d) I f+ G Lq{R^), f_,fdiv G Lq{R^),d€W^iR^), 
/div G g = (gi,... ,gjv+i) G h = {h,, ... G W^{R^)}, 

Xq = {F = (F+i, F_ 1 , F_2, F_3, F_4, Fi, F2, F3, F4, F5, Fe) \ F±i G Lq(R^), 
F_2,F_3,F_4 GL,(R(!), Fi,F2,F3,F4,F5 GL,(R'^), Fe€W^(R^)}. 

Then, there exist a constant Aq > 0 and operator families 

A±{X) €Rol{i:,,Xo,C{Xq,W^{R^))), iP_ gHo1(E,,Ao,^(^„VF,'(K-))), 
H(A) gHo1(E,,Ao,£(^„%"(K'^))) 
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such that u± = ^±(A)Fa, 7 r_ = 7 ^_(A)Fa and H = 7^(A)Fa are unique solutions of 
problem en) and (ESI) for any A e i;e,Ao and F = (f+, f_,/dw , fdiv , g, h, d) e Xg, 
and we have 

'^C(Xq,L^{m.f_)){{{'^^'rfG\A±{\) I A e Se^Aol) < C (£ = 0, 1), 
7^£(Ar,.L,(B-))({(Ta.)'V7^_(A) I A e E,.a„}) < c {1 = 0, 1), 
T^C(X„,W^(S.f^)){{{'^^TYG\%{\) I A e Se^Aol) < C (£ = 0, 1) 
with some constant c. Here, 

GiAl±(A) = (AAl±(A),Ai/ 2 vAl±(A),VM±(A)), GlUiX) = {m{\) ,Vn{\)), 

Fa = (f+, f_, A^/Vdiv , V/div , Afdiv , Ai/2g, Vg, Ah, A^/^vh, V^h, d), 

= {tt. g L,doc(M'^) I V^_ G Lg{R^)}, 

Hol(C/, X) denotes the set of all holomorphic functions defined on U with their 
values in X, V = {di,..., On) and = {didj \ i,j = 1,..., N). 

Remark 5. (1) F±i, F_ 2 , F_ 3 , F_ 4 , Fi, F 2 , F 3 , F 4 , F 5 and Fq are correspond¬ 
ing variables to f±, A^/^/dw , V/dw , fdiv, Ai/2g, Vg, Ah, Ai/^Vh, V^h and d, 
respecitively. 

(2) The proof of Theorem 15.31 is given in [T^ . 
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